Chromosome Oscillations in Mitosis 
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Successful cell division requires a tight reg- 
ulation of chromosome motion via the activity 
of molecular motors. Many of the key players 
at the origin of the forces generating the move- 
ment have been identified, but their spatial and 
temporal organization remains elusive[lj. The 
protein complex Kinetochore on the chromosome 
associates with microtubules emanating from one 
of the spindle poles and drives the chromosome 
toward the pole0, HI- Chromokinesin motors 
on the chromosome arms also interact with mi- 
crotubules, ejecting the chromosome away from 
the pole[^, In animal cells, a mono-oriented 
chromosome (associated to a single pole) period- 
ically switches between phases of poleward and 
away from the pole movement [6, |3J |£g, a behavior 
tentatively explained so far by the existence 
of a complex switchin g mech anism within the 
kinetochore itself[7, 9, 10, ll|. Here we show 
that the interplay between the morphology of 
the mitotic spindle and the collective kinetics 
of chromokinesins can account for the highly 
non-linear periodic chromosome motion. Our 
analysis provides a natural explanation for the 
origin of chromosome directional instability and 
for the mechanism by which chromosomes feel 
their position in space. 

The characterization of the forces leading to chro- 
mosome motion is one of the central questions in mi- 
tosis research The oscillatory movement of mono- 
oriented chromosomes observed in prometaphase that 
persists during chromosome congression, metaphase and 
early anaphase^ 0, H @i is a signature of the forces 
acting on the chromosome. The poleward (P) force ex- 
erted by the kinetochore is thought to be due to cytoplas- 
mic dyneins[2| and microtubule (MT) depolymerization 
in the kinetochore Q ■ Chromokinesin motors associated 
with the chromosome arms|l2j move toward the plus end 
of MTs @ and generate the force moving the chromosome 
away- from-the-pole (AP)[l[l3|. All the available models 
for the oscillatory movement of the chromosome (often 
referred to as "kinetochore directional instability" ||) 
have as a common feature that the kinetochore some- 



how controls the switching between P and AP phases 
JUS] However, the collective behavior of molec- 
ular motors can give rise to dynamical instabilities 
which have been observed in biological and biomimetic 
systems |l5l llij . Indeed, inhibiting chromokinesin ac- 
tivity leads to the disappearance of the oscillations and 
to the collapse of mono-oriented chromosomes onto the 
centrosome Q. In order to precisely assess the role of 
chromokinesins, we analyze the balance of forces on a 
mono-oriented chromosome (Fig. la). 




FIG. 1: Schematic representation of a mono-oriented chromo- 
some, (a) Sketch of a MT aster (dark green) interacting with 
a single chromosome (light green). The kinetochore (red) is 
connected to the pole through a bundle of MTs. Chromoki- 
nesin motors on the chromosome arms may be bound to a MT 
(dark blue dots) or unbound (light blue dots) . (b) Forces driv- 
ing chromosome motion: the kinetochore poleward force Fk 
(red), the polar ejection force created by the bound chromoki- 
nesins Fap (dark blue) and the friction force opposing motion 
— £r (black). The chromosome postion relative to the pole (or- 
ange arrow) is r and its velocity is r. (c) Binding/unbinding 
kinetics of chromokinesin motors, with rates kb and k u re- 
spectively. In the bound state, chromokinesins move toward 
the plus end of MTs with a velocity V. 
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The chromosome motion occurs at length and veloc- 
ity scales for which inertial effects are negligible (low 
Reynolds number). As a result, a difference between P 
and AP forces induces a viscous motion of the chromo- 
some, characterized by a phenomenological friction coef- 
ficient £. Force balance then reads 

Fap -F K -£f = (1) 

where r = dr/dt is the chromosome velocity (t being 
the time), Fk is the P force applied by the kinetochore, 
and Fap the AP force (Fig. lb). The friction coeffi- 
cient £ generically characterizes the resistance to chro- 
mosome motion, and pretends to encompass such effects 
as hydrodynamic friction, MT polymerization and de- 
polymerization and the activity of molecular motors in 
the kinetochore [Til llif . Chromosome micromanipulation 
during anaphase in grasshopper oocytes shows that the 
chromosome velocity f is independent of chromosome po- 
sition and decreases line arly with an external force oppos- 
ing the poleward motion[l9( ■ These observations strongly 
suggest that both the kinetochore force Fk and the chro- 
mosome friction £ are nearly constant. We adopt such 
minimalist description of the kinetochore activity here 
and we estimate F K ~ 700pTV and £ ~ 6 1(T 2 N.s/m 
from 

It has been shown that the kinetochore does not push 
on the chromosome, for instance through MT polymer- 
ization [jflEll ]. Here we consider that the AP force is cre- 
ated solely by the binding and displacement of chromoki- 
nesin motors on the MT aster (Fig. la). Chromokinesins 
stochastically attach to and detach from MTs with aver- 
age binding and unbinding rates kb and k u respectively. 
At a given time, only some amount n(t) of the total num- 
ber TV of chromokinesins permanently attached to the 
chromosome arms is bound to MTs and able to partici- 
pate to the AP force. Generically, the time evolution of 
the number n of bound motors may be written as 

^ = h(N - n) - k u n (2) 
dt 

The higher the concentration of MTs, the easier it is 
for a motor to find an attachment site, hence the bind- 
ing rate kb is larger in dense regions of the MT aster. 
In a monopolar spindle (Fig. la), the MT concentration 
decreases away from the pole, so the binding rate kb is 
a decreasing function of the chromosome position r. It 
is important to note that neither the P or the AP forces 
depend explicitly on chromosome position. All spatial 
information is contained in the binding rate fcfc(r) and 
reflects the morphology of the MT spindle. 

Once bound, both the motor velocity V and its un- 
binding rate are strongly influenced by the motor load. 
If the n bound chromokinesins are independent from one 
another, they equally share the total ejection force Fap 
and each motor feels a load Fap/ti. The rate of mo- 
tor unbinding increases exponentially with an applied 



load k u = ku^ exp[FApa/nkBT] (Kramers theory [22]|h 
where k^ is the unbinding rate at vanishing load, a is 
a phenomenological length and ksT is the thermal en- 
ergy. This exponential sensitivity to applied force has 
been observed for conventional kinesin [23j (for which 
k^ ~ 0.5s[24| and a ~ 1.3nm[2^) and also for myosin 
motors [2^. The velocity of a motor decreases with 
a force opposing motor movement, and vanishes at a 
particular stall force f s . For conventional kinesin, the 
force-velocity relationship is nearly linear, with a veloc- 
ity at vanishing force Vq ~ 0.6fim/s and a stall force 
f s ~ 6pjV|26l |27| . We adopt the linear relationship 
V = Vb(l — Fap /nf s ), with no substantial influence on 
our results. 

Identifying the chromosome velocity r with the 
chromokinesins velocity V on MTs and combining the 
equations above, we obtain a self-contained set of two 
coupled equations for {r, n} (see Supplementary Notes). 
The origin of chromosome movement is fairly clear. Close 
to the pole, the MT density is high, many chromoki- 
nesins attach and produce a large force that moves the 
chromosome away from the pole. Far from the pole, 
MTs are scarce, chromokinesins detach and the chro- 
mosome moves poleward due to the kinetochore force. 
Somewhere in between, there exists a fixed point where 
the systems remains still. It corresponds to a number 
of bound chromokinesins n s = Fk/ f s that exactly bal- 
ances the kinetochore force and to a chromosome position 
r s (given by kb(r s ) = k^e^n s / (N — n s )) where chromoki- 
nesin attachment and detachment fluxes exactly compen- 
sate (where / = f s a/kE,T quantifies the influence of the 
motor load on its detachment rate). 

The fixed point may be stable, in which case the chro- 
mosome remains at a fixed position r s , or unstable, lead- 
ing to permanent chromosome oscillations (see Supple- 
mentary Notes). Linear perturbation analysis around 
{r s ,n s } shows that the dynamical state of the chromo- 
some is controlled by three parameters. The dimension- 
less factor / = fga/ksT and the numbers n s = Fk/ f s 
and 7i£ = £Vo// s normalized by the total number of 
chromokinesins N. The fixed point is actually unsta- 
ble for a wide range of parameters satisfying < 
fn s (l — n s /N) — n s (Fig. 2). Note that the boundary 
between stable and unstable regimes does not depend on 
the precise MT distribution in the aster, provided the 
density decreases away from the pole (see Supplementary 
Notes). A reduction of the total number N of available 
chromokinesins has the effect of moving the system along 
a straight line in the parameter space (arrow in Fig. 2), 
and eventually to exit the oscillatory regime. The dis- 
appearance of chromosome oscillations has indeed been 
observed upon a drastic reduction of chromokinesins Q. 
We predict that the oscillation should stop for a precise 
number of chromokinesins: 2V C = fn 2 s /((f — l)n s — n^). 

Numerically computed chromosome motion in the un- 
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Stable states 




FIG. 2: Dynamical regimes of chromosome motion. Param- 
eter space representing the region of stable chromosome po- 
sition and chromosome oscillations (/ = 2). Increasing n s 
and correspond to increasing kinetochore force and chro- 
mosome friction, respectively. Decreasing the total chromoki- 
nesin number N (arrow) eventually leads to the disappearance 
of oscillations at N = N c (filled circle) . 



stable regime (Fig. 3a) display the characteristic sawtooth 
shaped oscillations observed in-vivo Q (Fig. 3c), indicat- 
ing that the system switches suddenly between phases 
of constant velocities. Indeed, the period of the oscilla- 
tion (~ min) is controlled by the viscous motion of the 
chromosome, while the switching between phases occurs 
over much shorter time scales, characteristic of motor 
binding/unbinding (<~ sec). Close to the pole, the high 
density of MTs results in the binding of a large number of 
motors (Fig. 3b), driving the chromosome AP at a veloc- 
ity close to their maximum velocity Vq. Thus, we argue 
that the velocity of the AP motion is a direct quantita- 
tive estimate of the chromokinesin velocity at vanishing 
load (Vb ~ 2/im/min in Newt lung cells pHl])- As the 
chromosome moves away from the pole, the density of 
MTs decreases and eventually reaches a value at which 
the attachment flux is too low to compensate the motor 
detachment. The remaining motors then detach rapidly 
(Fig. 3b) and the chromosome switches to P movement. 
The P phase occurs with almost no motors attached and 
the chromosome moves toward the pole with a constant 
velocity —Fk/£- The cycle is completed when the chro- 
mosome reaches a region of high enough MT density, 
where many motors abruptly attach and eject the chro- 
mosome. The ratio of AP and P velocities, which charac- 
terizes the symmetry of the oscillations, is approximately 
given by n^/n s . In Newt Lung cells (Fig.3c)Q, the pole- 
ward velocity is of order 2/xm/min, and the oscillations 
are roughly symmetrical. From their amplitude and pe- 
riod, we estimate N ~ 1000 and N c ~ 600. Hence, the 
inactivation of about half the chromokinesins would be 
sufficient in this case to suppress oscillations. 

We propose that chromokinesin binding onto MTs al- 
lows the chromosome to feel its position in space via 
the aster-like morphology of the MT spindle and that 




FIG. 3: Chromosome oscillations. Position of the chromo- 
some (a), and number of bound chromokinesins (b) as a func- 
tion of time. The theoretical analysis reproduces and explains 
the sawtooth shape of in-vivo oscillations ((c) from |T^). The 
chromosome oscillates between phases of poleward (P) and 
away-from-the-pole (AP) motion with nearly constant veloc- 
ities. The sudden switch from P to AP and AP to P phases 
corresponds to a dramatic binding (B) and unbinding (U) 
of chromokinesins respectively. The parameters are chosen 
to reproduce the amplitude and period of the in-vivo oscilla- 
tions in (c). For the sake of simplicity, the MT distribution 
in the aster is assumed to be isotropic, so that k^{r) ~ 1/r 2 . 
The AP velocity is an estimate of the chromokinesin velocity 
at vanishing load Vo = 2/j,m/min (see text). The remain- 
ing parameters are k^ = Is -1 , kb(r) = 30/(r(/xm)) 2 s _1 ), 
n s /N = 0.115, n ( /N = 0.052 and / = 2. 

chromokinesin unbinding from MTs is the force-sensitive 
mechanism at the origin of the chromosome instability. 
Within our framework, the highly non-linear, sawtooth 
shaped oscillations result from the combination of these 
two effects, and not from a mechanism inside the kine- 
tochore. The kinetochore might also be force-sensitive, 
e.g. via the activity of cytoplasmic dvneins|2^|. Nev- 
ertheless, spatial information is essential for chromosome 
oscillations. We argue that it is provided by the position- 
dependent chromokinesin attachment rate and the spa- 
tial information provided by the MT density in the aster. 

Finally, our analysis has implications for the congres- 
sion of bi-oriented chromosomes 13J. The symmetry of 
a bipolar spindle insures that chromosomes correctly lo- 
cate at the mitotic plate if the net force (Fk — Fap) 
toward each pole increases with the distance to the pole. 
Ostergreen|2Jj proposed that the kinetochore force Fk 
increases linearly with the length of the MT fiber con- 
necting it to the pole. We argue that the increase of the 
net poleward force is indeed nearly linear, but is rather 
due to decreasing polar ejection forces Fap away from 
the center of the MT aster (see Supplementary Discus- 
sion) . In conclusion, chromosome positioning may be 
due solely to the aster-like microtubule distribution and 
the kinetics of chromokinesins. 

We thank J. Prost for critical reading of the manuscript 
and J-F. Joanny, F. Jiilicher and J. Prost for stimulating 
remarks. O.C. thanks the European Network PHYNECS 
(HPRN-CT-2006-00312) and Ministerio de Education, 
Cultura y Deporte for financial support. 
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Supplementary Notes 

Dynamical regimes of chromosome movement 

Mono-oriented chromosomes in a microtubule (MT) 
aster are subjected to poleward (P) force exerted by one 
of the chromosome kinetochores and to an away-from- 
the-pole (AP) force produced by chromokinesin motors 
on the chromosome arms. The unbalance between these 
two forces leads to a viscous motion of the chromosome 
with a velocity r (r is the distance from the pole). The 
number n of chromokinesins exerting the AP force at 
one given time depends on the kinetics of motor attach- 
ment/detachment to MTs, itself dependent of the mor- 
phological properties of the MT aster. The chromosome 
movement is fully determined by a set of coupled dynam- 
ical equations for n and r, and the temporal evolution of 
the system corresponds to a given trajectory in the {n, r} 
"state space". Using the equations given in the main 
text, the set of coupled differential equations (dynamical 
system) reads: 



h = k b (r)(N -n) - k^ 



exp 



n s + nt 

/— — - I n 

n + 



r = V n 



n t 



(SN-1) 



All spatial information is contained in the binding rate 
kb(r), which includes the MT density in the aster. The 
force-sensitive mechanism leading to the instability arises 
from the force-dependent detachment rate and the col- 
lective dynamics of motors: k u = k^ exp^FApa/nksT). 
The ratio of the kinetochore force to the chromokinesin 
stall force n s = Fk / f s is the number of bound motors 
necessary to compensate the kinetochore force, and the 
number = £Vq/ f s is the ratio of the chromosome fric- 
tion to an effective chromokinesin friction. 

This dynamical system has a single fixed point defined 
by n = and r = 0, and given implicitly by {n = n Sl r = 



5 } where 



k b (r s ) = k° u e f 



N-n s 



(SN-2) 



The linearized dynamics for the perturbations {Sn, Sr} 
near the fixed point reads: 



d_ 

dt 



5n 
Sr 



N 



Vo 
n s -\-7lfi 



N-r 



dk b {r s ) 
Or 




(SN-3) 



A 



The matrix A specifies the linearized dynamics of the 
system around the fixed point. The stability of the fixed 
point is obtained from the eigenvalues of A, or equiva- 
lently from its trace and determinant: 



Det(A) 



■Vo 



Tr(A) = k°e f 



N - n. 



dk b (r s ) N -n s 
dr n s + n£ 
n, 

— 1 



+ / 



(SN-4) 



The determinant is always positive if enough motors are 
available (N > n s ) and for a decreasing concentration 
of MTs away from the pole (d r k b (r) < 0). The stability 
of the fixed point is fully specified by the trace of A. If 
Tr(A) < 0, the fixed point is stable and any perturbation 
relaxes back to the fixed point. If Tr(A) > 0, pertur- 
bations are enhanced, driving the system away from the 
fixed point. The equation Tr(A) = sets the boundary 
between stable and unstable regimes. The stability of 
the systems is entirely determined by the three param- 
eters: n s /N, n^/N and /. The chromosome is unstable 
for < n| = fn s (l — n s /N) ~ n s . The parameter / 
quantifies the influence of the load on the motor detach- 
ment rate. If / < 1, n| is always negative, implying that 



the fixed point is always stable. In other words, there 
exists a motor sensitivity threshold to force below which 
the chromosome always reach a stable position. 

In Fig lSN-fl we represent the lines along which h = 
(n— nullcline) and r = (r— nullcline). The chromoki- 
nesin kinetics is typically much faster than the chromo- 
some motion. As a consequence, the number of bound 
motor quickly adapts to a given chromosome position and 
the system always tries to follow the n- nullcline (n = 0). 
Fig lSN-lk shows a typical trajectory toward a stable 
fixed point. When the n- nullcline is non-monotonous 
fFig lSN-Jb l. the region of positive slope of the n-nullcine 
is unstable because the detachment flux of bound motors 
is too high to be equilibrated by the attachment flux. The 
linear stability analysis only specifies the transition be- 
tween stable and unstable regimes. Fie lSN-lb shows that 
once the system becomes unstable, the trajectories in the 
"state space" {n, r} evolve toward a limit cycle of the 
non-linear system, corresponding to permanent periodic 
chromosome oscillations. This can be shown formally us- 
ing the Poincare-Bendixon theorem (see J. Guckcnhcimcr 
& P. Holmes Nonlinear Oscillations, Dynamical Systems, 
and Bifurcations of Vector Fields (Applied Mathematical 
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FIG. SN-1: Stable and oscillatory trajectories in the "State 
space" of the dynamical system. Each point corresponds to 
a given number of bound motors n and chromosome position 
r (kfr (r), equivalently) . The lines of vanishing chromosome 
velocity (f = - dashed line) and of balance between motor 
binding and unbinding (n = - black line) intersect at the 
fixed point (filled circle). The phase trajectories (thick gray 
line) start from an (arbitrary) initial condition (open circle), 
(a) Stable regime: the n-nullcline is monotonous, and the sys- 
tems evolves toward the (stable) fixed point, (b) Oscillatory 
regime: the fixed point is unstable and the system follows 
a limit cycle corresponding to periodic oscillations between 
the extreme positions r< and r> . Starting near the pole, the 
phases of the oscillatory chromosome movement are: (1) AP 
motion with many bound chromokinesins, (2) sudden motor 
unbinding far from the pole, (3) P motion with no bound 
motors and (4) sudden motor binding near the pole. Parame- 
ters: / = 2, n s = 0.115AT and (for (a)) n e = 0.152iV (for (b)) 
n ( = 0.052iV. 



Sciences Vol. 42), Springer- Verlag, 1990). As can be seen 
in Fig lSJN-fb . the unstable range for the number of bound 
motors corresponds to a range of chromosome positions 
(r< < r < r>) between which the chromosome oscillates. 



Supplementary Discussion 

"Micromanipulation" of mono and bi-oriented 
chromosomes 

Molecular forces responsible for chromosome motion 
in mitosis may be probed by micromanipulation experi- 
ments where an external force is applied to the chromo- 
some. Here we investigate the effect of such force on the 
position of the chromosome within the mitotic spindle. 

The framework developed in the main text rely on the 
assumption of a constant force Fk exerted by the kine- 
tochore. If an additional external force F ext is applied 
to a mono-oriented chromosome in the direction of the 
pole, it merely shifts the effective value of the kinetochore 
force to Fk + F ext and modifies the chromosome motion 
accordingly. Bi-oriented chromosomes experience com- 



peting kinetochore forces from each pole and away-from- 
the-pole (AP) forces from both microtubule (MT) asters. 
The kinetochore forces being constant in our framework, 
they cancel out for each chromosome position. On the 
other hand, polar ejection forces implicitely depend on 
the chromosome position in the spindle via the chromoki- 
nesin binding rate and the local MT density. AP forces 
cancel at the mitotic plate due to the bipolar symmetry 
of the mitotic spindle (Fig lSD-lk ). Note that bi-oriented 
chromosomes may also display the oscillatory behavior 
described in the main text. Here, we assume for the sake 
of clarity that a bi-oriented chromosome reaches a sta- 
ble position, e.g. due to an increased friction from the 
bi-polar attachment. Chromokinesins on stable chromo- 
somes always exert their stall force f s and the spatial de- 
pendence of the AP forces arises only from the variation 
of the number of bound motors. A chromosome under an 
external force F ext directed toward the left pole is shifted 
leftward to a position r s satisfying the force balance: 



[ne(r s 



(r s ))fs 



(SD-1) 



where n r and rig are the number of chromokinesins bound 
on MTs from the left and right asters, respectively. 

In a monopolar spindle, the spatial variation of the 
chromokinesin binding rate follows the spatial distribu- 
tion of microtubules (MTs) in the aster and is given by 
the function fcb(r). In a bipolar spindle with a distance 
L between the spindle poles, we define the dimension- 
less function of local MT density a{r) = kf,(r)e~^ 
where r is the distance from the left pole. At the sta- 
ble position r s , equilibration of motor attachment and 
detachment fluxes on each aster leads to: 



ng(r) = n r (L — r) = N 



a(r) 



1 + a(r) + a(L — r) 



(SD-2) 



where N is the total number of available chromokinesins. 
Combining the two equations above gives a relationship 
between the applied force and the stable chromosome 
position: 



F,, 



a(r s ) - a(L - r s ) 

1 + a(r s ) + a(L - r.) ~ Nf s 



(SD-3) 



which could be directly measured in micromanipulation 
experiments. 

Eq l|SD-3|) is represented on Fig lSD-lb for various 
aster-like MT distributions. At vanishing external force, 
the stable position is at the mitotic plate r s = L/2. Due 
to the symmetry of the bipolar spindle, the chromosome 
displacement is linearly proportional to the external force 
even for fairly large displacements from the mid-plane. 
Mathematically, this arises from the existence of an in- 
flexion point for the left-hand-side of Eq (|SD-3J| . so that 
the quadratic term in the Taylor expansion near the mid- 
plane vanishes. Fig lSD-lb shows that the linear relation- 
ship is valid almost up to values of the external force 
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FIG. SD-1: (a) Number of chromokinesin motors bound to 
microtubules from the left (pink) and right (turquoise) asters. 
The balance of polar ejection forces and external force (black) 
is also represented for two symmetrical chromosome positions, 
(b) Stable chromosome position as a function of the exter- 
nal force. The chromosome displacement with respect to the 
spindle mid-plane varies linearly with the force almost up to 
chromosome collapse onto one of the poles. Curve are for pow- 
erlaw decreases of microtubule density with distance from the 
pole a(r) = r _/3 with j3 — 1.5 (red) f3 = 2 (green) and (3 — 3 
(blue) 



at which the chromosome collapses onto one centrosome: 
F ext ~ Nf s . According to our estimate, this force is in 
the nN range. 



The situation we described has to some extent been re- 
alized experimentally in (Hays, T.S., Wise, D. & Salmon 
E.D. Traction force on a kinetochore at metaphase acts 
as a linear function of kinetochore fiber length J. Cell 
Bio. 93 374-382 (1982)) by analyzing the metaphase 
stable position of multivalent chromosomes having more 
than two kinetochores. In this case, the excess number 
of kinetochores generate an extra force toward one of the 
poles. The resulting chromosome displacement with re- 
spect to the mid-plane has indeed been reported to be 
proportional to the excess number of kinetochores. 



